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ABSTRACT 



We consider the production of gravitinos and moduli fields from quantum vacuum fluctuations 
induced by the presence of scalar metric perturbations at the end of inflation. We obtain 
the corresponding occupation numbers, up to first order in perturbation theory, in terms of the 
power spectrum of the metric perturbations. We compute the limits imposed by nucleosynthesis 
on the spectral index n s for different models with constant n s . The results show that, in certain 
cases, such limits can be as strong as n s < 1.12, which is more stringent than those coming 
from primordial black hole production. 
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1 Introduction 



In the framework of the inflationary cosmology, density perturbations are generated when metric 
and inflaton quantum fluctuations become super-Hubble sized during inflation and reenter the 
horizon as classical fluctuations in the radiation or matter dominated eras jlj]. Typically, the 
scale corresponding to the present horizon size {Hq 1 ~ 3000 Mpc) left the horizon about 60 
e-folds before the end of inflation, and therefore one expects a primordial spectrum of metric 
perturbation spanning a range of scales from thousands of megaparsecs down to e~ m /H ~ 
10 -23 Mpc. The current observations of cosmic microwave background anisotropics and large 
scale structure have started providing us with precise measurements of the primordial spectrum 
at large scales (1000 - 1 Mpc). Smaller scales are much more difficult to probe since their 
evolution has already become non-linear or the corresponding fluctuations have been erased 
by various damping effects in the early universe. However, still there is a window open to 
those scales which is provided by primordial black holes. When a perturbation crosses inside 
the horizon with relative large amplitude 1/3 < 5 < 1, the perturbed region stops expanding 
and starts collapsing eventually forming a black hole || 0, ||. If such black holes do not 
evaporate fast enough they could overdose the universe and, therefore, limits on their primordial 
abundance are required. If they do evaporate by the present time, limits are still needed in 
order for the evaporation products do not disrupt primordial nucleosynthesis. Such limits can 
be translated into bounds on the power spectrum, which are particularly useful to constrain 
models of inflation in which perturbations grow at smaller scales, i.e., when the corresponding 
spectral index is larger than that of Harrison and Zeldovich n s — 1. 

Apart from the generation of primordial black holes, metric perturbations at small scales 
can have other interesting effects, which we will explore in this paper. The presence of metric 
inhomogeneities is known to break the conformal flatness of the Friedmann- Robert son- Walker 
(FRW) metric. In fact, such effect is more important as we move to smaller and smaller scales, 
and thus, in principle, one expects to find a strong deviation from conformal flatness right 
at the end of inflation, provided the smallest scales reenter the horizon with sufficiently large 
amplitude. However, we will find that even perturbations with very small amplitude can have 
important cosmological effects. One of the physical consequences of conformal non-invariance 
is the possibility of creating particles from vacuum fluctuations [f| |]. Recently 0, § this 
mechanism has been applied to the generation of large scale magnetic fields after inflation. 
Also in ||, the production of fermions from metric perturbations during preheating has been 
considered. Other sources for fermion production after inflation have been studied in JTO] . In 



this paper we will be interested in the creation of a different kind of relics which can also have 
important consequences in cosmology, namely gravitinos and moduli fields. 

In supergravity theories, the gravitino is the spin-3/2 superpartner of the graviton field. 
Its couplings to the rest of matter fields are typically supressed by the Planck mass scale, 
which implies that gravitinos can live very long and even decay after nucleosynthesis destroying 
the nuclei created in this period. A similar effect is due to the presence of moduli fields, 
corresponding to the radii of extra dimensions in higher dimensional theories. This imposes very 
strong contraints on the primordial abundance of this kind of particles, thus typically we have 
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n/s < 1(T 12 - 10~ 14 for masses of the relics in the range m = 10 3 - 10 2 GeV |JTj, [□], 0, 0, ^ • 
Traditionally, these constraints are translated into upper limits on the reheating temperature 
of the universe. In order to avoid their thermal overproduction we need Tr < 10 8 — 10 9 
GeV. More recently, non-thermal production of gravitinos and moduli during preheating has 



been considered [[L6[ [17], |TI| [19], ^TJ. In those works, it has been shown how the coherent 
oscillations of the inflaton field after inflation could produce a large amount of relics that could 
conflict in some cases with the above limits. In fact, in the case of moduli fields, the limits 
imposed on the reheating temperature can be as strong as T R < 100 GeV, if the gravitational 
or the production during preheating are taken into account | 20|| . In this work, we will show 



how in certain models of inflation, metric perturbations can also give rise to a relevant amount 
of gravitinos and moduli. In those cases, the nucleosynthesis limits discussed above impose 
stringent constraints on the perturbations power spectrum. 

The plan of the paper goes as follows. In Section 2, we consider the gravitino equations of 
motion in an inhomogeneous background and obtain their perturbative solutions. In Section 
3, using the Bogolyubov technique, we compute the total number of gravitinos produced as 
a function of the power spectrum and compare the results with the nucleosynthesis bounds. 
Section 4 is devoted to a similar analysis but using moduli fields. Section 5 contains the main 
conclusions of the work and finally we have summarized some useful formulae in the Appendix. 



2 Gravitino field equations in inhomogeneous backgrounds 

The massless Rarita-Schwinger equation is conformally invariant, which implies that gravitinos 
are not produced from vacuum fluctuations in a FRW background. In the case of massive grav- 
itinos in FRW, they can be produced either gravitationally |23| or due to the time dependence 
of the mass during preheating [T1], [I7|, [Tj| . However, since in this work we are interested in the 
effect of metric perturbations, we will consider that the mass of the gravitino is constant and 
very small and effectively we can neglect it in the calculations. This is a good approximation for 
models in which the gravitino mass is smaller than the Hubble parameter after inflation, since 
gravitinos will be produced with energies much higher than their mass. For simplicity, we will 
also consider only the production of helicity ±3/2 particles which is technically simpler. The 
helicity ±1/2 production exhibits additional subtelties related to the fact that their equations 
of motion are only consistent if the gravitational background is a solution of the supergravity 
equations of motion. In order to get consistent inflationary supergravity models, one usually 
requires the introduction of several chiral superfields. In that case, the equations of motion for 
the helicity ±1/2 gravitino cannot be decoupled from those of the fermionic components of the 
superfields ^T|. In addition, in the context of the production during preheating, if the infla- 



tionary sector is decoupled from the sector responsible for the present supersymmetry breaking, 
it is the fermionic partner of inflaton what is produced more abundantly rather than the lon- 



gitudinal gravitino ||22|| . These points make the calculations more involved. However, in our 
case, we do not expect any new physical effects coming from helicity ±1/2, since on physical 
grounds, the production mechanism is purely gravitational and the helicity ±1/2 components 
will be produced with similar abundance to the helicity ±3/2 states. 
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Let us then consider the massless Rarita-Schwinger equation in an external gravitational 
background: 

e^WvD^ = 0. (1) 

where ip a is a Majorana spinor satisfying ip a = CipJ with C = 27 2 7° the charge conjugation 
matrix. The covariant derivative is given by: 

D^ = {d, + Y^ a -Yl^ x (2) 

and the spin connection by: T M = — |r° fe [7 a , %]. With these definitions, we have D^ u = and 
because of the totally antisymmetric tensor, the Christoflel symbols do not contribute to (U). 
As usual, latin indices a,b, . . . refer to the tangent-space tensors, whereas greek ones fi, is, . . . are 
used for curved background objects. The different gamma matrices are related by 7 M = e^7 a , 
where the vierbein satisfies: e a ^e b y g^ u = r] ab , with rj a b the Minkowski space metric. 

For the background metric we will take the following form: 

9iw = 9% + V ( 3 ) 

where 

g^dx^dx" = a 2 (r])(dr] 2 - <%dxW) (4) 

is the flat FRW metric in conformal time and 

h^dx^dx" = 2&a 2 (ri)(dr] 2 + 5ijdx l dx j ) (5) 

is the most general form of the linearized scalar metric perturbation in the longitudinal gauge 
and where it has been assumed that the spatial part of the energy-momentum tensor is diagonal, 
as indeed happens in the inflationary or perfect fluid cosmologies [0. In this expression $(77, x) 
is the gauge invariant gravitational potential. 

Contracting equation ([[]) with 7a 7 M we get: 

2iD A (7 CT Vv)-2iWA = (6) 

As commented before we are asumming that metric perturbations are classical perturbation 
which are produced during inflation. In this sense, it is a good approximation to consider that $ 
vanishes asymptotically in the past. Since we are interested in the smallest scale perturbations, 
which are going to produce the leading effects, we will consider only those perturbations that 
reenter the horizon during the radiation era. In this case, once they reenter they start oscillating 
with damped amplitude, for that reason, we will also take $ — > when rj — > 00. The vanishing 
of the perturbations allows us to define conformal vacuum states in the asymptotic regions. 
Thus, in those regions, the above equation reduces to a Dirac-like equation % fl^ip^ = (see 
[[Uj) where the index (0) denotes the unperturbed object. The corresponding positive frequency 
solution with momentum p and helicity A = ±3/2 can be written as: 
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j u(p,±)e,(p,±)e"* (7) 



l2pVa{rj) 

where the polarization vectors are given by 

t^(p,+) = -^(0, COS6 1 cos0 — 2 sin 0, cos6* sin ± i cos0, — sin#) (8) 
v2 

€fj,(p,—) = /=(0, cos 6* cos0 + i sin 0, cos 6* sin — i cos0, — sin#) (9) 

with p^ = (p, p sin 9 cos 0, p sin 9 sin 0, p cos 9) and the normalization e*(p, m)e M (p, n) = —5 mn , 
p^e^ip, m) = p^e* (p, m) = 0. The helicity r, s = ±1/2 spinors are chosen such that they satisfy: 
w(p,r)u(p, s) = 2p5 rs . Here we are working in a finite box with comoving volume V and we 
will take the infinite volume limit at the end of the calculations. Notice that the above solutions 



satisfy the additional constraints |JT6 j 



p^f^ix) = 0, yvi 0,wt (ar) = 0. (10) 



We will look for perturbative solutions of equations (^) in the form: 

^ = < 0)P ' A + ^ 1) + ... (11) 

where, as mentioned above, ip^ p ' x is the solution of the unperturbed equation given in ([7p. 
Expanding the equations @ up to first order in the perturbation we find: 

2Wf{ 1 ^ !1 ^) - 2% ^ (0 Vi 1} - 2i ^ (1 Vi 0) = (12) 

where we have used the constraint equations (p~Q|) . In Appendix A we have given the expressions 
for the perturbative expansions of the different terms in the above equation. Since we are only 
interested in the evolution of those states which asymptotically (in the past and in the future) 
give rise to helicity ±3/2 gravitinos, we will project these equations along those helicity states. 
The projectors in the asymptotic regions are given in Fourier space by 

P£ 3/2 (k)=e»*(k,±)P ±1/2 (k) (13) 

where P r are helicity r = ±1/2 projectors satisfying P r (k) u(k,s) = u(k,s)S rs . Let us assume 
for simplicity that the three momentum is pointing along the z-direction k^ = (\k\, 0, 0, A;). In 
this case, the projectors take a very simple form: 

e"*(fc,±) = i=(0, 1,^,0) 

P ±1/2 (£) = i(l± 7 5 7°7 3 ) (14) 
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Notice that the \i = equation in ( |T2"D does not contribute to the projected equations. Using 
the fact that P^i 2 {k)k^ip = and P^^ik)^^ = for an arbitrary spinor tp, we see that the 
first term in ( |T2"D vanishes when projected and we are left with: 



PLs/2(k) 



d 3 x 
(2vr) 3 / 2 



(15) 



Using the formulae in Appendix A we get the following expression for the equations of 
motion: 



where 



^% 2 (k, V )=a 1 / 2 P^ 3/2 (k)4 1 \k,v) 



(16) 



(17) 



with ip^ (k, T]) the corresponding Fourier mode, defined as usual by f(k, 77) = (2ty) 3 / 2 / d 3 xe lkx f(x) . 
We have already defined ?/>f^ A (p, r/) in (|7|). Finally we have denoted <£>(£ + p, rj) simply by $. 

Notice that since (p, 77) is a solution of the unperturbed equations of motion, the sec- 
ond term in the previous equation vanishes. Finally, in order to reduce (|Tj]) to a standard 
(inhomogeneous) harmonic oscillator equation, we multiply by (7°<9o — irfkj), and we get: 



(«9 2 + A; 2 ) M) + J ±3/2 

The spinor current is given by: 
1 







J. 



±3/2 



(p2 + k 2 - 2^ 7 °7 J > - *7°(%7 J - Pl°)& P±3/2(M° )X (P, V) 



where prime denotes derivative with respect to conformal time rj. Because of the presence of the 
inhomogeneous current, the initial positive frequency solution, with momentum p and helicity 
A in (^) will evolve into a linear superposition of positive and negative frequency modes, with 
different helicities and different momenta. Thus in the asymptotic future we find: 



^\x) 



Tj — >00 



J2 12 



A'=± k 



u(k, A / )e /J (fc, A') ^(kg-kn) 
2kVa(r]) 



+ 



u^k,X%(k,X') 

Ppkxx' i — e 

l2kVa(r]) 



kri) 



(19) 
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In order to obtain the Bogolyubov coefficients, we need to solve (|T8p. Up to first order in 
perturbations, we have: 



fe A /2 (M) =ft±l%\k,v) + l I" J4*i2*HK*-rf))drf (20) 



'/» 



where 7/0 denotes the initial time in the remote past when the perturbations were switched 
off. Comparing this expression with the above expansion, it is very simple to get an explicit 
expression for the Bogolyubov coefficients [p4 ]: 



P pkXT = -—— / Ti^(*,T)J±8/2C-*' (21) 



'/o 



where 971 denotes the present time or any other instant in time after inflation in which the 
perturbations vanish again. Notice that the helicity projectors P± 3 / 2 project positive frequency 
modes on ±3/2 helicity states, whereas negative frequency modes are projected on =f3/2 states, 
that is the reason why the (3 coefficients subindex are changed with respect to those of the 
current J±3/2- The total number of gravitinos with comoving momentum k is given by: 

N k = EE l/W'l 2 ( 22 ) 

XX' V 



3 Gravitino production from metric perturbations 

As commented before, we will concentrate only in the effect of super- Hubble scalar perturbations 
whose evolution is relatively simple. For single-field inflationary scenarios we have Jl|: 

^(k,Tj) = C k -4- (~ [ a 2 dri)+D k -^ (23) 



adrj \a 

the second term decreases during inflation and can soon be neglected. Thus, it will be useful to 
rewrite the perturbation as: <§>{k,rj) = CkJ-(r)). During inflation or preheating, super-Hubble 
perturbations evolve in time, whereas they are practically constant during radiation or matter 
eras. Notice that this expression is not appropriate in two-field models of inflation, since in that 
case super-Hubble gravitational fluctuations can be exponentially amplified during preheating 



2^] due to parametric resonance. 

The power spectrum corresponding to (1231) is given by: 



= ^ = * (£) s ~ (24) 

For simplicity we have assumed a power-law behaviour with spectral index n s and we have set 
the normalization at the COBE scale \q — 3000 Mpc with As — 5 ■ 10~ 5 . Although there are 
some models which predict this kind of behaviour, in general we will have a dependence of the 
spectral index on the scale. However, in order to compare our results with those obtained from 
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black holes, we will keep this simple form. In principle, there will be a maximum frequency cut- 
off k which corresponds to the perturbation with the smallest size produced during inflation. 
Typically this scale is roughly determined by the size of the horizon at the end of inflation 
k ~ ajHj where the I index denotes the end of inflation (for a detailed discussion see M). 
We can obtain an explicit expression for the occupation number (|22|) in terms of the power 
spectrum. Taking the continuum limit Yl p ~^ (2n)~ 3 / 2 V J d 3 p, we get: 



where we have used kn <C 1 which is valid for super-Hubble modes and 6 is the angle between 
p and the z axis. The total comoving number density of gravitinos is given by N = / d 3 kN^, 
typically this integral is dominated by the upper limit of integration, so that we can approximate 
iV ~ k 3 N k . Notice again that we are assuming no enhancement of high- momentum modes 
during preheating as it happens for instance in two-field models of inflation. It is then enough 
to calculate N K : 



Because of Pauli exclusion principle < 1, the violation of this bound signals the breakdown 
of the perturbative approach. 

In order to get the ratio n/s we need to calculate the entropy density at the end of inflation. 
If we assume that all the energy density in the inflaton field is instantaneously converted into 
radiation, we get a very high reheating temperature, which is already excluded by the thermal 
production of gravitinos as commented before. Therefore, for realistic models, reheating and 
thermalization should occur sufficiently late, so that we can have low reheating temperature. 
In that case, we must consider the stage of inflaton oscillations at the end of inflation. If the 
inflaton potential close to the minimum behaves as V ~ <p a , then the energy density during 
oscillations scales as p ~ a -Ga/(a+2) ||26|| . Thus, although the physical number density scales as 
n ~ a~ 3 , the entropy density will do as s ~ p 3 / 4 ~ a _9a /( 2 ( a + 2 )). This implies that during the 
oscillations and depending on the value of a, the number density n can be diluted with respect 
to the entropy density, it could remain constant or even increase. In the following, we will 
consider the usual potentials with a = 2,4 (notice that we are approximating the potential by 
simple powers only close to the minimum, during inflation their behaviour can be completely 
different). For a = 2, the equation of state of the oscillating scalar field can be considered as 
that of non-relativistic particles, in this case n/s decreases in time. For a = 4, the oscillations 
behave as radiation and n/s is constant. Thus, taking this into account and setting k = ajHj, 
we get for the ratio at the end of reheating: 



XX' J 



(2tt) 3 / 2 



dp dVt 





(26) 




(27) 



where M P = (8nG) 1//2 . In order to compute N K in the previous equation we need to know the 
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Figure 1: Abundance of gravitinos n/s versus spectral index for a model with a = 4 and 
Hi = 10 13 GeV (full line), Hj = 10 10 GeV (dashed line) and Hj = 10 7 GeV (dotted line). The 
results are independent of the reheating temperature 



value of the ratio n/kc- If the Hubble parameter remains constant during inflation we have: 
ft/kc = CLi/a>c, where a>c denotes the scale factor at the moment when the scale Ac left the 
horizon. Therefore, in terms of the number of e-folds of inflation, we can write: 



!L _ p N(k c ) 

k c 



In the case a = 2, N(kc) is given by the well-known expression P7J : 

N(k c ) = 53 + - In Tr + - In ■ 



3 10 w GeV 3 10 u GeV 



(28) 



(29) 



where V is the value of the inflaton potential during inflation. A similar expression can be 
obtained for a = 4: 



N(k c ) = 56 + In 



v 1 / 4 



W u GeV 



(30) 



which is independent of the reheating temperature. Substituting back in ( P7Q and using 
and (EBD, we get: 



n V2^A% ( Hj 



3/2 



4 — a 
2a 



y/ZHiMp, 



exp(K - l)N{k c )) 



(31) 



In Fig.l we have plotted n/s as a function of the spectral index for different models with 
Hj = 10 13 , 10 10 and 10 7 GeV with a = 4. In this case the results are independent of the 
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Figure 2: Limits on the spectral index ft $ clS db function of the gravitino mass 7713/2, from the 
generation of massive unstable gravitinos 
duction of D and 4 He 
D and 3 He O 



The right curve corresponds to the hadropro- 
The left curve to the photofission of D and the photoproduction of 
The model of inflation has a = 4 and Hj = 10 13 GeV 



reheating temperature. The nucleosynthesis bounds on the spectral index coming from the 
effects of the hadronic and radiative decays of gravitinos on elemental abundances is summarized 
in Fig. 2 for a typical model with Hj = 10 13 GeV. We see that the strongest bound comes from 
gravitino masses around 100 GeV for which n s < 1.12. This limit improves that imposed by 
primordial black holes production. In fact, for Tr < 10 9 GeV, such limit comes from black 
holes which are evaporating today, and is given by n s < 1.28 |3|, |J. Notice that for n s = 1.12 
and n/kc ~ 10 26 (which is the range of scales spanned in the model considered) we have 
V<&(k) — 3 • 10~ 6 <C 1, i.e. we are well inside the perturbative region. Notice that such 
perturbations are unable to create black holes, but still we have an important effect coming 
from the creation of relics. 

In Fig. 3 we plot the number density versus spectral index for a model with a = 2 and 
Hi = 10 13 GeV, for different values of the reheating temperature. Now because of the relative 
growth of entropy during the inflaton oscillations, the production is much weaker than in the 
previous case. In fact, in order to get a cosmologically relevant abundance with low reheating 
temperature, the spectral index should be larger than 1.4. However for those values we find 
Vq> > 1, i.e. we are out of the perturbative regime. For that reason for a = 2 we cannot obtain 
limits on the spectral index. In Fig. 3 we have also plotted the relative abundance for a model 
with high reheating temperature (Tr = 10 15 GeV) with dashed line. The results show that 
even if we ignored the thermal production of gravitinos, we will find again the problem, but 
now coming from metric perurbations. 

In [0 a weak production was obtained with numerical calculations for the V = m 2 <j) 2 /2 and 
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Figure 3: Abundance of gravitinos n 3 / 2 /s versus spectral index n s for a model with a = 2 and 
reheating temperature Tr = 10 9 GeV (full line), Tr = 10 15 GeV (dashed line) and Tr = 10 5 
GeV (dotted line), for R I = 10 13 GeV 

V = A0 4 /4 models of chaotic inflation. This is due to the fact that those models predict a 
power spectra with negative tilt, i.e. n s < 1. 



4 Moduli production 



String theory and other higher-dimensional models include scalar fields that parametrize the 
shape of extra dimensions and whose couplings to matter fields are Planck mass suppressed. 
This kind of fields can give rise to different cosmological problems ||29|| . First, if they were 
displaced from their minima in the early universe, we find the so called classical moduli problem 
in which the energy density of the oscillations of the scalar fields behaves as non-relativistic 
matter, and can give relevant contribution to the energy density of the universe. If these fields 
decay during nucleosynthesis, we find a problem similar to that produced by gravitinos. The 
typical bounds on the primordial abundance of moduli are also similar n/s < 10 -12 — 10~ 14 
PH . Second, moduli field can be produced from quantum fluctuations during inflation. In this 
case, the breaking of conformal invariance comes either from the minimal coupling to the scalar 
curvature or from their own mass. Thus the typical equation of motion for moduli is given by: 



□ + ml + iR) X = 



(32) 



In general, the moduli mass acquires corrections, coming from supersymmetry breaking effects 
ion or from non-renormalizable effects in 
where C is a model dependent parameter. 



during inflation or from non-renormalizable effects in the superpotential, thus one has m 2 x 



m 



2 + C 2 H\ 
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Figure 4: Abundance of moduli n x /s versus spectral index n s for a model with a = 4 and 
Hi = 10 13 GeV (full line), if/ = 10 10 GeV (dashed line) and Hj = 10 7 GeV (dotted line). The 
results are independent of the reheating temperature. 



It has been shown |30| that the classical moduli problem is not present in the cases C <^ 1 
or C ^> 1. The gravitational quantum production of moduli could also be avoided if conformal 
invariance is recovered, i.e., ifC = 0, £ = 1/6 and m = 0, although such fine tuning in the 



parameters seems to be difficult to obtain without invoking additional symmetries [|1J|. Since 
we are interested in the effect of metric perturbations on the moduli production, we will consider 
in the following that they are the only source of conformal symmetry breaking. In that case 
the above equation of motion reduces to: 



(33) 



Using the form of the metric in @ and the formulae in the Appendix, we find the following 
form of the linearized equation of motion: 



x" - didoc - 4$'x' - ^(mm - $"x = o 

where x = a X- Again we look for solutions in the form x — + X^ + • • • where 

1 



x (0) 



D i(px-pri) 



(34) 



(35) 



corresponding to the initial conformal vacuum state. Fourier transforming the equation, we 
find for the first order correction: 



^ i y , (k, v ) + k 2 x^(k, v ) + j = o 



(36) 
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Figure 5: Abundance of moduli n x /s versus spectral index n s for a model with a = 2 and 
reheating temperature Tr = 10 9 GeV (full line), Tr = 10 15 GeV (dashed line) and Tr = 10 5 
GeV (dotted line), for E I = 10 13 GeV. 



where 



(37) 



Using the analogous result to the gravitino case, we find for the Bogolyubov coefficients: 



Pi 



pk 



and finally for the occupation number we get: 

d 3 p \C p+k \ 2 



r Je- ikri dr] 

J r?n 



(38) 



(2tt) 3 / 2 AkpV 



AipF' + -(p + k) 2 F- T" ) dq 

vo v 3 



(39) 



Assuming again that the total number density is dominated by those modes with k ~ k, we 
have: 



327T° f p° f 
Nk - fo Y3/2 / dP^^K) \ 



(27t) 3 / 2 



IT' + ^— T - ^— ) drj 
V 12p \,») 



(27T) 



3/2 



72 



(40) 



where the time integral has been estimated to be ~ n/(12p), since T ~ O(l)- 



The ratio n/s can be calculated in a similar fashion using (|2~7 ) and the expression we have 
just obtained for N K . The results are plotted in Figs. 4,5 for different models with a = 2,4. 
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In the case in which there is no relative dilution of n with respect to s, we obtain that for 
a moderate value of the nucleosynthesis bound n/s < 10~ 13 , the spectral index should be 
n s < 1.18, whereas for the strongest one n/s < 10 -14 , we get n s < 1.14 for Hi ~ 10 13 GeV. 
These limits are slightly weaker than those coming from gravitinos, in part this is due to the 
absence of spin degrees of freedom in this case. For a = 2 the limits are much less stringent as 
expected. 



5 Conclusions 



In this paper we have studied the production of helicity ±3/2 gravitinos and moduli fields 
from metric perturbations after inflation. The results show that those metric perturbations 
with very small wavelengths reentering the horizon right at the end of inflation induce a strong 
deviation from conformal flatness which is responsible for the production of a non-negligible 
amount of relics. In particular, if the power spectrum has a positive tilt n s > 1, as predicted 



by some models of hybrid inflation ]3T| , the production could conflict with the limits imposed 
by nucleosynthesis. This in turn allows us to set stringent constraints on the spectral index 
in the particular case in which n s does not depend on the scale. When the energy density in 
the inflaton oscillations scales as radiation, the limits that we obtain n s < 1.12 are stronger 
than those coming from primordial black hole production. In the case in which the energy in 
oscillations scales as non-relativistic matter, the bounds are much weaker. When the spectral 
index depends on the scale, it is possible to translate the nucleosynthesis bounds into limits on 
the power spectrum at the smallest scale Using the results in the previous sections, we 

obtain < 10~ 6 , in models of inflation with Hi ~ 10 13 GeV and a = 4. In conclusion, 

we consider that the production of gravitinos and moduli fields provides us with a new tool 
which can be useful to study the region of very small wavelengths in the primordial spectrum 
of metric peturbations. 
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A Useful formulae 



In this Appendix we show some useful results for the perturbative expansions of different 
geometrical objects, (see also ||). 
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Vierbein 

The vierbein expansion corresponding to the metric tensor in (||) is given 
where 



h° = 2$, h) = -2$<5V 



Gamma matrices 

The curved gamma matrices can be expanded as: 



where 



7m = el% = 7 i 0) + 7? } + 



Spin connection 

For the spin connection appearing in the fermionic derivatives we find: 

r M = rj» + r a) + ... 

p(O) _ « r (0) _ (0) (0) 

h 8 V p p a ^ a 

= -^„[7f,7 ( °»1+47f»7f 
Christoffel symbols 

The non-vanishing Christoffel symbols for the metric (0) are given by: 



pA 


_ r (0)A , r (l)A , 

(IV ~ [IV ~ ' ' 




F (0)0 
1 00 


a> r (o)o a' 
a J a 




F (0)i 
1 Oj 






r (i)o 

1 00 


* ) 1 0« 




r (i)o 




r (l)i 

1 oj — 


r (l)< 






r (i)< 

1 00 
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Scalar curvature 

Finally, for the scalar curvature we find up to first order: 

R = — [6a" - 12$ a" - 24a' $' + 2a^^$ - 6$" a] 

(47) 

Spinors 

The normalized spinors with helicity ±1/2 and momentum p that we have used in the text 



are: 



e -^/2 cos (|) \ f _ e -«A/2 sin ^ 



(4* 



e i<t>/2 gm 



in (|) j ' y e^cos(f) 
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